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A MULTILEVEL CORRECTION SCHEME FOR NONSYMMETRIC 
EIGENVALUE PROBLEMS BY FINITE ELEMENT METHODS 


HEHU XIE* AND ZHIMIN ZHANGt 

Abstract. A multilevel correction scheme is proposed to solve defective and nodefective of 
nonsymmetric partial differential operators by the finite element method. The method includes 
multi correction steps in a sequence of finite element spaces. In each correction step, we only need 
to solve two source problems on a finer finite element space and two eigenvalue problems on the 
coarsest finite element space. The accuracy of the eigenpair approximation is improved after each 
correction step. This correction scheme improves overall efficiency of the finite element method in 
solving nonsymmetric eigenvalue problems. 
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1. Introduction. As we know, the numerical approximation of eigenvalue prob¬ 
lems plays a central role in the analysis of the stability for nonlinear partial differential 
equations. For example in fluid mechanics, the analysis of the hydrodynamic stability 
always leads to a nonsymmetric eigenvalue problems (see [l][7l|8]). The stability of the 
underlying flow depends on the real part of the eigenvalue which has the smallest real 
part (see [DIE]). For more details, please refer [T1[5J[TH]. The aim of understanding 
the stability of nonlinear partial differential equations naturally leads to the compu¬ 
tation of the eigenvalue problems with some numerical methods. The main content 
of this paper is to design an efficient finite element method to compute nonsymmetric 
eigenvalue problems. 

Recently, a multigrid method is designed to solve the self-adjoint eigenvalue prob¬ 
lem based on a type of multilevel correction method [iniiiiiiiiiiH]- But as we know, 
the analysis of the stability for nonlinear partial differential equations always leads to 
nonsymmetric eigenvalue problems mm and the extensions of the multilevel method 
for self-adjoint eigenvalue problems to the nonsymmetric ones is not direct [21 HSl US] 
and needs more analysis. So the purpose of this paper is to propose a multilevel cor¬ 
rection scheme to solve nonsymmetric eigenvalue problems based on the finite element 
method. In the past, a two-grid finite element method was proposed and analyzed by 
Xu and Zhou in |23j for symmetric eigenvalue problems. Latter, Kolman used this idea 
to design a two-level method for nonsymmetric eigenvalue problems in [2] . Yang and 
Fan |25j also studied a two-grid method for nonsymmetric eigenvalue problems. As 
an alternative approach, in niiiiiii], the authors used a recovery technique PPR 
to improve the convergence rate for both symmetric and nonsymmetric eigenvalue 
problems. All these methods are designed for the nonsymmetric eigenvalue problems 
under the assumption that the ascent of the concerned eigenvalues is only one which 
means the algebraic eigenspace is the same as the geometric eigenspace. 

Along the line of multilevel correction method, here we present a multilevel correc¬ 
tion scheme to solve nonsymmetric eigenvalue problems without the ascent assump¬ 
tion. With the proposed method, solving nonsymmetric eigenvalue problems will not 
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be much more expensive than solving corresponding source problems. The correction 
method for eigenvalue problems in this paper is based on a series of finite element 
spaces with different levels of accuracy which are related to the multilevel method (cf. 

ED]). 

The standard Galerkin finite element method for nonsymmetric eigenvalue prob¬ 
lems has been extensively investigated, e.g. Babuska and Osborn E1I3], Chatelin 

and references cited therein. Here we adopt some basic results in these papers 
to carry on error estimates for our multilevel correction scheme. It will be shown 
that the convergence rate of the eigenpair approximations can be improved after each 
correction step. 

Our multilevel correction procedure can be described as follows: (1) solve an 
eigenvalue problem in the coarsest finite element space; (2) solve a source problem in 
an augmented space with the associated eigenfunction from (1) as the load vector; (3) 
solve the eigenvalue problem again on a finite element space constructed by enhancing 
the coarsest finite element space with the eigenfunction obtained in step (2). Then 
go to step (2) for the next loop. 

In this method, we replace solving the eigenvalue problem in finer finite element 
spaces by solving a series of boundary value problems in a series of nested finite element 
spaces (with the finest space as the last one) and a series of eigenvalue problems in 
the coarsest finite element space; and yet, we achieve the same accuracy as solving 
the eigenvalue problem in the finest space. It is well known that there exist multigrid 
methods that solve boundary value problems with the optimal computational work 
(cf. |21jb Therefore, combined with the multigrid method, our correction method 
improves overall efficiency in solving nonsymmetric eigenvalue problems (cf. |181I19] 1. 

An outline of the paper goes as follows. In Section 2, we introduce the finite ele¬ 
ment method for nonsymmetric eigenvalue problems. An one level correction scheme 
is described and analyzed in Section 3. In Section 4, we propose and analyze a mul¬ 
tilevel correction algorithm to solve nonsymmetric eigenvalue problems by the finite 
element method. Some numerical examples are presented in Section 5 to validate our 
theoretical analysis and some concluding remarks are given in the last section. 

2. Discretization by finite element method. In this section, we introduce 
some notation and error estimates of the finite element approximation for nonsym¬ 
metric eigenvalue problems. Throughout this paper, the letter C (with or without 
subscripts) denotes a generic positive constant which may be different at different oc¬ 
currences. For convenience, we use symbols <, >, and «, such that xi < yi,X 2 ^ 2 / 2 , 
and xs ~ ys have meanings: xi < Cij/i, X 2 > C 22 / 2 , and C 3 X 3 < 2/3 < C 3 X 3 , for some 
constants Ci,C 2 ,C 3 , and C 3 that are independent of mesh sizes (cf. [20]). 

We consider the following eigenvalue problem: 

Find A S C and u such that 

{ —V • (AVu)-I-b • Vu-I-()>it = X(pu, in H, 

u = 0, on on, (2.1) 

J^ip\u\^dn = 1 , 

where H C Ti‘^ is a bounded polygonal domain, A G b € (j) is a function 

defined on and (/? is a real positive function with (p > cq >0. 

We define V := Hq{H) with the usual norm || • ||i. The corresponding variational 
form of m can be stated as follows: 

Find {X,u) gC xV such that 6 (m, u) = 1 and 

a{u,v) = Xb{u,v), \fv G V, 


( 2 . 2 ) 
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where 


a{u, v) = {AVu, Vv) + (b • Vm, v) + v), 

b(u, v) = {(pu, v) 

with (•, •) denoting the inner product in the space The corresponding adjoint 

eigenvalue problem is: 

Find (A,m*) gC xV such that b(u*,u*) = 1 and 

a{v,u*) = Xb{v,u*), Wv € V. (2-3) 

In the sequel, we also use the norm ||?;||{, = \/b(v, v) which is equivalent to the L^{n) 
II 0 - Here the bilinear form a(-, •) is assumed to satisfy 

a{w, v) 


norm 


^ sup 


v&V Fill 


and ||ui||i < sup f. ’ ^ , Vic G V. 
vev Fill 


(2.4) 


We further assume that a(-,-) is V-elliptic, i.e., 

||m||^ < Re a(M, u), \/uGV. (2-5) 

2.1. Operator reformulation. We introduce the operators T, T* G C{V) de- 
hned by the equation 

a{Tu,v) = b{u,v) = a{u,TtfV), \/u,vGV. (2-6) 

The eigenvalue problem (12.211 can be written as an operator form for A ^ 0 (denoting 
p := A-i): 


Tu = pu, 

(2.7) 

T^:U* = pu* 

(2.8) 


for the adjoint eigenvalue problem (12.31) . Note that ellipticity condition (12.51) guaran¬ 
tees that every eigenvalue A is nonzero. It is well known that the operators T and 
T* are compact. Thus the spectral theory for compact operators gives us a complete 
characterization of the eigenvalue problem (12.21) . 

There is a countable set of eigenvalues of (12.211 . Let A be an eigenvalue of problem 
(1^ . There exists a smallest integer a which are called the ascent such that 

7V((r-M)“)=iV((T-^)“+i), (2.9) 

where N denotes the null space and we use the notation /i = A“^. Let M{\) = 
Mx.^ = and Q{X) = Qa./x = N{T—p) denote the algebraic and geometric 

eigenspaces, respectively. The subspaces Q{X) C M(A) are finite dimensional. The 
numbers m = dimM(A) and q = diin(5(A) are called the algebraic and the geometric 
multiplicities of p (and A). The vectors in M(A) are generalized eigenvectors. The 
order of a generalized eigenvector is the smallest integer j such that (T — pYu = 0 
(vectors in Q{X) being generalized eigenvectors of order I). Let us point out that a 
generalized eigenvector of order j satishes 


a{uYv) = Xb{uYv) + Xa{u^ ^,v), 'iv GV, 


(2.10) 
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where ^ is a generalized eigenvector of order j — 1. 

Similarly we define the spaces of (generalized) eigenvectors for the adjoint problem 


= = and Q*{\) = = N{T, - fl). 

Note that fi is an eigenvalue of T (A is an eigenvalue of problem (12.21) 1 if and only 
if /2 is an eigenvalue of T* (A is an eigenvalue of adjoint problem (I2.3|) l with the ascent 
a and the algebraic multiplicity m for both eigenvalues being the same. 

2.2. Galerkin discretization. Now, let us define the finite element approxi¬ 
mations for the problem (12.21) . First we generate a shape-regular decomposition of 
the computing domain fl C TZ‘^ {d = 2, 3) into triangles or rectangles for d = 2 (tetra¬ 
hedrons or hexahedrons for d = 3). The diameter of a cell iF € 77i is denoted by Hk- 
The mesh diameter h describes the maximum diameter of all cells K € Th- Based 
on the mesh Th-, we construct a finite element space denoted hy Vh V. In order 
to define our multilevel correction method, we start the process on an initial mesh 
I'm with mesh size H and the initial finite element space Vh defined on Tb. In this 
paper, the finite element space 14 is assumed to satisfy 


llw'/illi sup 

VheVh 


a{wh,Vh) 


j II II ^ aivh.Wh) 

and ||in/i||i 4 sup —^- 

vh(iVh 


'iwh e 14- 


( 2 . 11 ) 


The standard Galerkin discretization of the problem (12.21) is the following: 
Find {Xh,Uh) €C xVh such that b(uh,Uh) = 1 and 


a{uh,Vh) = Xhb{uh,Vh), Vv/j e 14- (2-12) 

Similarly, the discretization of the adjoint problem (12.31) can be defined as: 

Find (XhjU’l) GC xVh such that b{u*^,uD = 1 and 

a{vh,ul) = Xhb{vh,ul), yvh e 14- (2-13) 

By introducing Galerkin projections Ph, Ph ^ ^h) with the following equa¬ 

tions 

a{PhU,Vh) = a(u,Vh), Vu € V, Mvh G 14, 
a{vh,u) = a{vh,PY)^ G V, Mvh G 14, 

the equation (12.121) can be rewritten as an operator form with ■= X^ (Note that 
Ph is a bounded operator). 


PhTuh = fihUh- (2-14) 

Similarly for the adjoint problem (I2.13|) . we have 

= flhul (2.15) 


2.3. Spectral approximation of compact operators. Let ^ be an eigenvalue 
(with algebraic multiplicity m) of the compact operator T. If T is approximated by 
a sequence of compact operators Th converging to T in norm, i.e., lim ||T — T/i || i =0, 

h-s-O-l- 

then for h sufficiently small /i is approximated by m numerical eigenvalues 
(counted according to their algebraic multiplicities) of Th, i.e.. 


lim ^j^h = M for J = 1, • • • , ru- 

/t—)-0+ 
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The space of generalized eigenvectors of T is approximated by the subspace 

771 

Mu{\) = ^ V((T^ - (2.16) 

j=i 

where is the smallest integer such that 

We similarly define the space QhW = Qh^ = SjLi — fJ-j,h) and counterparts 
Qh{X) for the adjoint problem . 

Now, we describe a computational scheme to produce the algebraic eigenspace 
Mh{X) from the geometric eigenspace Qh{X) = • • • ,Uq^h} corresponding to 

eigenvalues • • • , Xq^h}, which converge to the same eigenvalue A. 

Starting from all eigenfunctions in the geometric eigenspace QhiX) (of order 1), 
we use the following recursive process to compute algebraic eigenspaces (cf. [16]) 

f a{u] h.Vh) - Xj^hKu^j^f^,Vh) = Xj^ha{u^~h,Vh), yvhGVh, , . 

1 b{u^^ f^,Vh) = 0 , yvhGQhiX), ^ ■ 

where £ > 2, ^ is the general eigenfunction of order £ and h — '^j,h € Qh{X) for 

j = 1, • • • ,9- 

With the above process, we generate the algebraic eigenspace 

corresponding to eigenvalues {Ai^^i, • • • , Xq^h, • ■ • , Xm,h}, which converge to the same 
eigenvalue A. Similarly, we can produce the adjoint algebraic eigenspace M^(X) from 
the geometric eigenspace Qi!j(A). 

For two linear spaces A and B, we denote 

0(yl,i?)= sup inf Ijw — u||i, ^{A,B)= sup inf ||w —'!;||b, 

weA,\\w\\i=i'"^B weA,\\w\\b=i'"^^ 


and define gaps between A and B in || • ||i as 


0(A, B) = max {0(A, B), &{B, A)}, 

and in || • ||b as 

(2.18) 

<P{A,B) = max{$(A,B),$(B, A)}. 

(2.19) 

Before introducing the convergence results of the finite element approximation for 
nonsymmetric eigenvalue problems, we define the following notation 

ShiX) = sup inf llu-u/^lli, 

ueM{\),\\u\\i=i'"heVh 

(2.20) 

Sl{X)= sup inf lilt*-u/, 111 , 

m*gm*(a).||«*||i=i 

(2.21) 

PhiX) = sup inf \\u-Vh\\b, 

ueMi\),\\u\\i=i'"heVh 

(2.22) 

pI{X) = sup inf ||u* - u/tllb. 

(2.23) 

Pa{h) = sup inf llTZ-ulli, 

/eFII/ll6=i'''^^'* 

(2.24) 
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VaW= sup inf ||r,/--!;||i. (2.25) 

In order to derive error bounds for eigenpair approximations in the weak norm || ■ ||{,, 
we need the following error estimates in the weak norm || • ||b of the finite element 
approximation. 

Lemma 2.1. m Lemma 3.3 and Lemma 3.4-]) 

77Q(h) = o(l), ? 7 *(h)=o(l) as0, (2.26) 

and 

Ph{X)<v:ih)SHi\), (2.27) 

PUX)<Vaih)SUX)- (2.28) 


The following theorem is a basic tool for our error estimates. 

Theorem 2.2. m Section 8]) When the mesh size h is small enough, we have 

0(M(A),M^(A)) < 5,,(A), 0(M*(A),M*(A)) < 5^(A), (2.29) 

$(M(A),M,(A))<p,(A), ^iM*iX),M*iX))<pUX), (2.30) 

\X-Xn\<SniX)SUX), (2.31) 

where A/i = ^ with Ai_h, • • • , Xm,h converging to X. 

3. One correction step. In this section, we present an one-step correction pro¬ 
cedure to improve the accuracy of the current eigenvalue and eigenfunction approxi¬ 
mations. This correction method contains solving some auxiliary source problems in 
a finer finite element space and two eigenvalue problems on a coarse finite element 
space. 

Assume that we have obtained the algebraic eigenpair approximations {Xj^hk j Uj^hk ) G 
TZ X Vh^ and the corresponding adjoint ones {Xj^h^ hk) £ ^ for j = i, • • • ,i + 

m — 1, where eigenvalues {Xj^hk}j~LT~^ converge to the desired eigenvalue Xi of (j2.2D . 
Now we introduce a correction step to improve the accuracy of the current eigenpair 
approximations. Let Vh^+i C 1^ be the conforming finite element space based on a 
finer mesh which is produced by refining Th^, in the regular way. We start from 

a conforming linear finite element space Vh on the coarsest mesh TA to design the 
following one correction step. 

Algorithm 3.1. One Correction Step 
1. For j = i, - ■ ■ ,i m — 1 Do 

Solve the following two boundary value problems: 

Find Uj^hk+i G such that 

= b{uj,h^,Vh^^,), (3.1) 

Find u* G '^hk+i such that 

I ^ '^'^hk+i £ ^hk+i- (3-2) 


End Do 
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2. Define two new finite element spaces 

VhM+i =Vh ® , • • • , Ut+rn-lM+l } 

and 

^HM+i =Vh ,wr+„^-i, 

Solve the following two eigenvalue problems: 

Find i\j,hk+i,Uj^hk+i) € 7?. X VH,hk+i such that b{uphk+i,Uj^hk+i) = 1 and 
a{uj,hk+i,VH,hk+J = >^phk+M'a]M+i^'^H,hk+i), '^VH,hk+i ^y^.h^+i- ( 3 - 3 ) 
Find e 7^ X such that , u* J = 1 and 

a{vHMk+i,ulh^^J = Xj,hk+ib{'aH,hk+i’'>^j,hk+i)’ '^vhm+i ^yHM+i- ( 3 - 4 ) 

3. Choose 2q eigenpairs {Xph^+i , Uj^hk+i YjtT^ and {Xj^hk+i, Yj^V^ de¬ 

fine two new geometric eigenspaces 

Qhk+I (-^i) = span{ui^/j;.^j ; ■ ■ ■ I Ui^q—I^hk+I } 

and 


=span{w*^,--- 

Based on these two geometric eigencpases, compute the corresponding alge¬ 
braic eigenspaces 


(Aj) = spanjMi,,,^^!, • • 

! Ui.^-^n—l,hk+l } 

(3.5) 

and 



^4+1 (A*) = spanjui,,!,^^,-- 

j Ui+m-l,hk+i }■ 

(3.6) 

The final output is: 



Correction(V h, {Xj^hk YjtT~^^^hk 

(A.),V4*JAO,14.+0. 



Remark 3.1. Since in Step 1 of Alaorithm 1 3. A the solving process for the bound¬ 
ary value problems is independent of each other for different j, we can implement them 
in parallel. Furthermore, the designing for this algorithm does not need the ascent as¬ 
sumption as in BW- 

Theorem 3.1. Assume there exist real numbers efn,{Xi) and ef,^{Xi) such that the 
given eigenpairs (A*), in One Correction Step \3.1\ have 

following error estimates 


0(M(AO,M,,(AO)<e/.,(AO, 

(3.7) 

e { M *{ X ,), M *^{ Xfi )< el ^{ X ,), 

(3.8) 

<i>(M(A.),M,,(A0)<r;:(i7)e,,(A0, 

(3.9) 
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<P{M*{X,),Ml{X,))<rj,{H)el^{X,). (3.10) 

Then after one correction step, the resultant eigenpair approximation 

have following error estimates 

0(M(A,), (A,)) < (A,), (3.11) 

0(M*(A,),M;:^^^(A.)) <4.+,(A.), (3.12) 

$(M(A0,A/4,^,(A0) < v:{H)eH,^,{X,), (3.13) 

3>(A/r(A.), A/4V^(A.)) < ^a(i?)4.+,(A.), (3.14) 


where Eh^^.iXi) ■= vl{H)£h^{X,)+Sh,,+,{X,) andel^^^{X,) := Va{H)£l^{X,)+Sl^^^{X,). 
Proof. From (12.101) . there exist the basis functions {“j}j=r ^ M{Xi) such that 

C z+m —1 \ 

Pifc(A^)Mfc,^' , Vu e F, (3.15) 

where Pjk{') denotes a polynomial of degree no more than a for k = i, - ■ ■ ,j with 
Z'li(Ai) = Ai and PjkiXi) = 0 for j<k<i + m — 1. We can define a matrix 
^ := [Pj+i — i)i<j— 1 e such that 

a{U,v) = b(VU,v), Vu € V, (3.16) 

where U := {ui, ■ ■ ■ ,Ui+m-i)'^■ It is easy to know that the matrix V is nonsingular 
providing Ai 0. 

For each from the definitions of 0(M(Ai), M/i^(Ai)) and $(M(Ai), Mh,^{Xi)), 

there exist a vector TZj := (ci, • • • , CmY' S such that 

\Wj,hk-T^^U\\i<ehk{\), for j = i, • • • + m - 1, (3.17) 

\\u 3 M-'^^U\\o<'ql{H)ehk{\), for j + m-1. (3.18) 

For any vh^+i G Vh^+i , we have 

- Ph,+Ji^V-^U,Vht,+f)\ = 

= h{uf - nJV-^VU,VH,^,) = |6(u^^'= - TZjU,Vk,^,)\ 

^ V*a(.H)£h^{Xi)\\vh,,+,\\i, foT j = i,--- ,i + m-l. (3.19) 

From ()2.1ip and (j3.19L the following estimate holds 

WUjm+i - Ph,^,njr-^uh < v:iH)ek,{x,), 

for j = i, • • • , i + m — 1. (3.20) 


Combining with the error estimate 

\\njv-^u - PH,^,nJv-^uh < 4.+,(4), 

for j = f, • • • , f + m — 1, (3.21) 


W'^jTk+i 


we have 


njr-^uh < iil{H)eh,{X,) + 4.+,(A.), 








Multilevel correction for nonsymmetric eigenvalue 


9 


ioi j = i, ■ ■ ■ ,i + m — 1. 

After Step 3, from the definition of Vn.hk+i p.22p . we derive 

sup inf \\u-VH,hi+A\i 

< sup inf \\u-Vh,+,\\i 

< sup inf will 

.-Tijr-^uh 


(3.22) 


< 


max 


j—i,-" ,i+m —1 


\\Uj,hk 


where := span{M'‘''+\ • • • , 


(3.23) 


Similarly, 


sup inf \\u*-VH,hk+A\i 


< 




(3.24) 


Then from the error estimate results stated in Theorem l2.2l for the eigenvalue problem 
(see, e.g., [31 Section 8]) and (I3.23I) - (I3.24I) . the following error estimates hold 


Q{M {Xi), Mhk^iiXi)) < riil{H)ehk[Xi) + 5hk+i{Xi), 
0(MyA),M^*^^^(A)) < r?,(iJ)e);jA)+<5^,^^(A). 

These are the desired estimates (13.111) and p.l2l) . Furthermore, 


(3.25) 

(3.26) 


<^{M{Xi),Mh^^,{Xi)) <rjl{H) sup inf \\u - vhm+iWi 

MeM(A).||M||l=l ’'H.fefc + l + l 




(3.27) 


where 


?j:iH):= sup inf ||T,/- < r;:(i7). (3.28) 

Then we obtain (|3.13|) . A similar argument leads to (|3.14|) . □ 

4. Multilevel correction scheme. In this section, we introduce a multilevel 
correction scheme based on the One Correction S'tep id.ll The method improves ac¬ 
curacy after each correction step, which is different from the two-grid methods in 

0113135]. 

Algorithm 4.1. Multilevel Correction Scheme 
1. Construct a coarse conforming finite element space Vh^ on Tki such that Vh C 
and solve the following two eigenvalue problems: 

Find {Xhi,Uhi) €TZx Vh^ such that h{uh^,Uh^) = 1 and 

a{uhi,Vh^) = XhM'U'hi,Vh^), Vu/ii G Vh^. (4.1) 

Find (XfinUl^^) G 77. x Vh^ such that = 1 and 

aivhi,Ki) = ^hib{vhi,ulj, yvh, G Vh^. 


(4.2) 
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Choose 2q eigenpairs {Aj M I Uj^hj Yjti ^ and {Xj^hi , YjYl ^ which approx¬ 
imate the desired eigenvalue Xi and its geometric eigenspaces of the eigen¬ 
value problem and its adjoint one Based on these two geometric 

eigenspace, we eompute the corresponding algebraic eigenspaces Mh^Xi) '■= 
space{ui,/ii, • • • , Ui+rn-iM } and {Xi) := space{u* , • • • , }. Then 

do the following correction steps. 

2. Construct a series of finer finite element spaces Vh^, • ■ • , on the seguenee 
of nested meshes Th^ i • j Thn (cf. mw- 

3. Do fc = 1, • • • , n — 1 

Obtain new eigenpair approximations {{Xj^hk+iYjtJT 
by Alaorithm \3.1\ 

Correction{VH, {A,.^ (A.), 

End Do 

Finally, we obtain eigenpair approximations {{Xj^hnYjYT~^^ ^h„iXi), M^^{Xi)). 

Theorem 4.1. After implementing Algorithm \4.1\ the resultant eigenpair ap¬ 
proximation {{Xj^hnYjYi'~^ T ^h„{Xi), M^^{Xi)) has following error estimates 


0(M(AO, Mh„ (A,)) < eh„ (A,), (4.3) 

4>(M(A0,M,„(A.)) < pl{H)ehftXY, (4.4) 

0(A/r(A.),M*„(A.))<4JA.), (4.5) 

$(M*(A.),M^*„(A.)) < pa{H)elS^i), (4.6) 

|A^,/^„-A,| <e/,„(A,)4„(A,), (4.7) 

where \h„ = ^ £h„{X^) = YYl=i 3*a{H)^~^5hk{XY and 

elSXY = Yl=iVa{H)-Hl^{XY. 

Proof First, the following estimates hold 

e{MiXY,MhftXi))<ehftXY, (4.8) 

<i>{M{XY,PhftXY) <v:{hi)euAXY <v:{H)euAXY, (4.9) 

0(M*(AO,M,;(AO) <4,(A^), (4.10) 

^M*{XY,M*^fXY)<Va{hi)eY{XY < Va{H)el{XY. (4.11) 

Then we set EhftXi) := ShftXi) and e))j(Ai) := A^^(Ai). 

By recursive relation and Theorem 13.11 we derive 

e{MiXY,MhftXY)<ehAXY=v:{H)eh„_ftXY+ShJX^) 


< v:{H)"eu„_AXY + v:{H)Sn„_AXY + ShftXA 

n 

(4.12) 

k^l 

and 

n 

<^{M{XA,MnftXA)<ii:{H)Y,pl{H)^-Hhk{XA. (4.13) 

fe=l 

These are the estimates (|4.3D and (14.41) and the estimates (14.51) and (14.61) can be proved 
similarly. From Theorem 12.21 (14.31) and (|4.5L we can obtain the estimate (14.7F □ 
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5. Numerical results. In this section, we give some numerical results to il¬ 
lustrate the efficiency of the multilevel correction scheme defined by Algorithm 14.11 
Here, we solve the following eigenvalue problem 

J — Au-|-b-VM = Xu, in H, 

\ u = 0 , ondn, > 

where b = G is a constant vector and H = (0,1) x (0,1). This example 

comes from [71|5]. We choose bi = 1 and &2 = 1/2 in Subsections 15.11 and 15.21 Then 
we choose bi = cos( 7 rxi) sin( 7 rx 2 ) and &2 = ~ sin( 7 rxi) cos( 7 rx 2 ) in Subsection 15.31 We 
also choose a complex vector b in the final example. 

When 6 i = 1 and 62 = 1/2, the problem (15.11) is nonself-adjoint, but all of its 
eigenvalues are nondefective (all algebraic eigenfunctions are of order 1 ) and real 
numbers 


Xk,e 


bl + bl 

4 


+ {e+f)TT\ 


(5.2) 


for k,(. G . 

The corresponding eigenfunctions can be chosen as real functions 

Uk^i = exp b 2 X 2 '^ sin(fc 7 rxi) sin(£ 7 ra; 2 ). (5.3) 

The corresponding adjoint eigenvalue problem has eigenvalues Xk,f: and eigenfunctions 

biXi -I- b2X2 \ 


= exp - 


-j sin(fc7ra;i) sin(.^7rx2). 


(5.4) 


5.1. Multi-space way. In this case, finer finite element spaces are constructed 
by increasing polynomial degrees of the beginning finite element space on the same 
mesh. We first solve the eigenvalue problem (j2.12ll by linear finite element on a rela¬ 
tively coarser mesh 7/r, then perform the first correction step with quadratic element, 
followed by cubic element for the second correction step and quartic element for the 
third correction step. Our initial mesh Th is obtained from the Delaunay triangula¬ 
tion followed by four levels of regular mesh refinement. Figure 15.11 depicts errors for 
the first eigenvalue (5/16-|-27r^) approximation, and Figure [521 plots numerical errors 
for the eigenfunction and the corresponding adjoint eigenfunction associated with the 
first eigenvalue. 

Furthermore, FigurejASlprovides numerical results for the summation of the errors 
for the first 6 eigenvalues: 5/16 -I- [27r^, Stt^, Stt^, Stt^, IOtt^, IOtt^]. 

From FiEfures l5. 1115.31 we find that each correction step improves the convergence 
order by two for eigenvalue approximation, and by one for eigenfunction approxima¬ 
tion when the exact eigenfunction is sufficiently smooth. 

To end this subsection, we make a comparison with the PPR method [13]. We 
see from Figure [ATI that the two-level correction scheme by the multi-space way has 
slightly better accuracy than the PPR method. However, the two-level correction 
needs to solve two extra boundary value problems while the PPR method only need 
to perform a local recovery at each node. Thus, we should say that the PPR method 
has better efficiency than the two-level correction under regular mesh refinement when 
the eigenfunction has regularity iJ^(f2) nlF^’°“(H). Nevertheless, three and four-level 
correction will outperform the PPR method. 








12 


Hehu Xie and Zhimin Zhang 



Fig. 5.1. Here, denote the eigenvalue approximation by linear element, is the eigenvalue 
approximation by the first correction with quadratic element, the eigenvalue approximation by 
the second correction with cubic element, A^ the eigenvalue approximation by the third correction 
with quartic element 




Fig. 5.2. Here, and denote the eigenfunction approximation and its adjoint approxima¬ 

tion by linear element, u^ and are eigenfunction approximation and its adjoint approximation 
by the first correction with quadratic element, and u^^, eigenfunction and its adjoint approx¬ 
imation by the second correction with cubic element, and u^^, eigenfunction and its adjoint 
approximation by the third correction with quartic element 


5.2. Multi-grid way. An alternative way of the multilevel correction scheme 
is to construct finer finite element spaces by mesh refinement. We first solve the 
eigenvalue problem (j2.12[l in the linear finite element space on an initial coarse mesh 
Th (Thi ■= Th)- Then we refine the mesh regularly with the resultant meshes Thk 
satisfying hk = 2^~^H for (A; = 2, • • • , n), and solve auxiliary source problems f|3.ip 
and (|3.2p in the linear finite element space Wfc defined on Th^ the corresponding 
eigenvalue problems (|3.3p and (|3.4p in VH,hk- We have the following estimate 

n n -| 

k=l k-1 

and similarly (A) ~ hn, which implies that the multilevel correction method 
achieves the optimal convergence rate if the initial mesh size H is reasonably small, 
say 77 = 1/4 as we will use in our numerical tests. 
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Fig. 5.3. Approximation errors for the summation of the errors for the first 6 eigenvalues 
by the multi-space way. Here, Aj ^ denote the eigenvalue approximation by linear element, Xj ^ 
is the eigenvalue approximation by the first correction with quadratic element, A| ^ the eigenvalue 
approximation by the second correction with cubic element, A| ^ the eigenvalue approximation by 
the third correction with quartic element 



Fig. 5.4. Comparison with the PPR method in m when bi = 10 and 62 = 1 . Xj^, eigenvalue 
approximation by linear element; X^, eigenvalue approximation by the first correction with quadratic 
element; Appr, the eigenvalue approximation by the PPR method 


Numerical results for the first eigenvalue A = 5/16 + 27r^ and the two associ¬ 
ated eigenfunctions are demonstrated in Figures 15.51 and 15.61 respectively. Here we 
use the uniform meshes with H = 1/4. Furthermore, Figure 15.71 provides nu¬ 
merical results for the summation of the errors for the first 6 eigenvalues: 5/16-1- 
pTT^, IOtt^, IOtt^] with H = 1/8 and H = 1/16, respectively. 

We observe from Figures ISTSlIS.71 that our multilevel correction method with the 
multi-grid way produces eigenvalue and eigenfunction approximations with the opti¬ 
mal convergence rate. Therefore, we can combine the multigrid method for boundary 
value problems and our multilevel correction scheme (cf. |111118] ) to achieve better 
efficiency for nonsymmetric eigenvalue problems. 

5.3. Eigenvalue problem on T-shape domain. In this subsection, we con¬ 
sider the eigenvalue problem (EH) on the L-shape domain H = (—1,1) x (—1,1)\[0,1) x 
(—1,0]. Since H has a reentrant corner, the singularity of eigenfunctions is expected. 
As a consequence, the convergence rate for the first eigenvalue approximation is 4/3 by 
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Fig. 5.5. Approximation errors for the first 
H = 1/4 



eigenvalue 5/16 + 2tt'^ by the multi-grid way with 



Fig. 5.6. Approximation errors for the first eigenfunction and its adjoint by the multi-grid 
way with H = 1/A 


the linear finite element method on quasi-uniform meshes. Since the exact eigenvalue 
is unknown, we choose an adequately accurate approximation A = 9.95240442893276 
as the exact first eigenvalue for our numerical tests. 

Our multilevel correction scheme is tested on a sequence of meshes Th {Thi ■= 
Th), Th2, • • • 1 Thr, produced by the adaptive refinement (cf. [13 El]). Here the ZZ 
recovery method (cf. [26] 1 is adopted as the a posteriori error estimator for eigenfunc¬ 
tion and adjoint eigenfunction approximations ^J\\uh — + ||it^ — ^*110 h- Figure 

15.81 shows the initial mesh and the mesh after 12 adaptive iterations. Figure [5^ gives 
the corresponding numerical results for the adaptive iterations. 

From Figure [5^ we observe that the multilevel correction method works well on 
adaptive meshes with the optimal convergence rate. Furthermore, the situation is 
very different from the two-gird [91E31E5] method in that the initial mesh has very 
little impact on the finest one. Thus the multilevel correction method can be coupled 
with the adaptive refinement naturally. 

5.4. Eigenvalue problem with complex vector. In this subsection, we test 
the multilevel correction scheme for the problem (ED with complex vector b = [1 -|- 
2i, 1/2 — i]"^. We use the multi-space and multi-grid ways as in Subsections 15.11 and 
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Fig. 5.7. Approximation errors for the error summation of the first 6 eigenvalues by the multi¬ 
grid way with 77 = 1/8 (left) and 1/16 (right) 


Initial mesh 


Mesh after 12 iterations 



Fig. 5.8. The initial mesh and the one after 12 adaptive iterations for the L-shape domain 


15.21 respectively, to check the multilevel correction scheme. Figure 15.101 shows the 
numerical results for the first 6 eigenvalues. It is observed from Figure [5.101 that the 
multilevel correction method defined in Algorithm 14 .1 1 can also work very well for the 
nonsymmetric eigenvalue problems with complex vector. 

6. Concluding remarks. In this paper, we propose and analyze a multilevel 
correction scheme to improve the efficiency of both defective and nondefective non¬ 
symmetric eigenpair approximations. In this multilevel correction, we only need to 
solve eigenvalue problems in the coarsest finite element space. Sometimes, we also 
need to compute the algebraic eigenspace based on the geometric eigenspace when 
the ascent is larger than 1. 

Furthermore, our multilevel correction scheme can be coupled with the multigrid 
method to construct a parallel method for eigenvalue problems (see, e.g, [TT] [121 
[T51 m]). It can also be combined with adaptive techniques (cf. [T7]) for singular 
eigenfunction cases. These will be our future work. 

A final remark. As long as higher eigenvalues are concerned, the multi-space 
way is preferred (than the multi-grid way). We can see it clearly by comparing 
numerical accuracies for summations of the first 6 eigenvalues in §5.1 and §5.2. 
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Fig. 5.9. Approximation errors of the first eigenvalue and the a posteriori errors of the 
associated eigenfunction and adjoint eigenfunction 



Fig. 5.10. Approximation errors for the summation of errors for the first 6 eigenvalues by the 
multi-space way (left) and the multi-grid way with H = 1/8 (right). Here, Aj ^ denote the eigenvalue 
approximation by linear element, A| ^ is the eigenvalue approximation by the first correction with 
quadratic element, ^ the eigenvalue approximation by the second correction with cubic element, 
A4 ^ the eigenvalue approximation by the third correction with quartic element 
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